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Abstract A new polarization survey of the northern sky at 1.41 GHz is presented. The observations were carried out using 
the 25.6 m telescope at the Dominion Radio Astrophysical Observatory in Canada, with an angular resolution of 36'. The 
data are corrected for ground radiation to obtain Stokes U and Q maps on a well-established intensity scale tied to absolute 
determinations of zero levels, containing emission structures of large angular extent, with an rms noise of 12 mK. Survey 
observations were carried out by drift scanning the sky between -29° and +90° declination. The fully sampled drift scans, 
observed in steps of 0.25° to ~ 2.5° in declination, result in a northern sky coverage of 41.7% of full Nyquist sampling. The 
survey surpasses by a factor of 200 the coverage, and by a factor of 5 the sensitivity, of the Leiden/Dwingeloo polarization 
survey I Sooelstra 1972) that was until now the most complete large-scale survey. The temperature scale is tied to the Elfelsberg 
scale. Absolute zero-temperature levels are taken from the Leiden/Dwingeloo survey after rescaling those data by the fac- 
tor of 0.94 . The paper describes the observations, data processing, and calibration steps. The data are publicly available at 
|http://www.m pifr-bonn.mpg.de/div/konti/26msurvey or http://www.drao.nrc.ca/26msurvey 
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1. Introduction 

The dete ction of linear polarization in the Galactic r adio 
emission dWesterhout et alJl962HWielebinski et al.ll962l) was 
the final proof that low-frequency radio emission is gen- 
erated in Galactic magnetic fie lds by the synchrotron pro- 
cess. Subsequent surveys (e.g. Berkhuii sen & Brouwl Il963t 
[wielebinski & Shakeshaft 1964) established the distribution 
of the polarized emission across the northern sky. The most 
prominent region of polarized emission, about 40° in extent, 
was found in the direction of / - 140° and b - 9° with no ap- 
parent corresponding count erpart in total iri t ensity . The south- 
ern sky was surveyed by Mathewson et alJ (1 19651) . These ob- 
servations with the Parkes telescope were made with a smaller 
beam (48') but the effective sampling was on a 2° grid. All 
the early surveys were made at the frequency of 408 MHz 
with rather low angular resolution, 2° in Dwingeloo and 7.5° 
in Cambridge. Over the following years polarization surveys 
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were extend e d to h igher frequencies, first to 610 MHz by 
iMuUer et ai] ( Il963h . later to 620 MHz by iMathewson etalJ 
( 1 19661) . and to 1407 MHz by[Bingham ( 196^! 

A major step forward in this field was the multifrequenc y 
mapping of the northern sky by Brouw & Spo elstral (Il976l) . 
Maps at 5 frequencies between 408 MHz and 141 1 MHz were 
presented. The absolute caUbration of the survey was related to 
the earlier observations of lWesterhout et alJ ( 1962) wi th a cor- 
rection of the polarization definition proposed by BerkhuiisenI 
(1975). However, at the highest frequency of 1411 MHz the 
sampling was not complete. The vectors were shown mostly 
on a 2° grid corresponding to the 408 MHz survey resolution. 

A major advance in polarization studies of the Galactic 
emission has been made through observations at higher angular 
resolution. The Effelsberg Medium Latitude Survey (EMLS) 
at 1.4 GHz (.Uvanikeret al. 1998. 1999; Reich et al.. 2004) has 
an angular resol ution of 9.4'. Th e Canadian Galactic Plane 
Survey (CGPS: 'Gray et al.' 'l998l Il999t iTavlor et alJ 120031: 
Uvaniker et al. 2003). and th e Southern Galactic Plane Survey 
(SGPS: Gaensl eret alJl200ll) have angular resolution of order 
1', also at 1.4 GHz. Valuable at these surveys are, all miss in- 
formation on the broadest structures. The CGPS and SGPS are 
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interferometer surveys, and complementary information on the 
broadest structure must be obtained with single-antenna po- 
larimetry; the CGPS undersamples structure larger than 45', 
and the SGPS structure larger than 30'. The EMLS is observed 
one rectangular area at a time, and the data reduction technique 
employed suppresses structures larger than about 5° in extent. 

Total intensity in the EMLS is restored to the correct zero 
level by extracting spatially filtered data from the Stockert sur- 
vey 1 llReichlll982tlReich & Reic h 1986). The same technique 
cannot be used for polarization data because the only avail- 
able dataset with an absolute zero lev el, the Leiden/Dw ingeloo 
polarization survey (LDPS hereafter; ISnoelstrJl97i . IS very 
sparsely sampled. There is a pressing need for fully sampled 
polarization data at 1 .4 GHz to complement the recent high- 
resolution mapping. 

A further motivation that deserves mention is the cur- 
rent interest in Galactic foregrounds to observations of the 
Cosmic Microwave B ackground (e.g. lBaccigalupi et alJl20oil 
iBuri gana & T , a Portal EnnZl Here it should be noted that 
Faraday effects, which play a minor role at frequencies of 
30 GHz and higher, are believed to affect most of the polarized 
emission at 1.4 GHz. Observed polarized intensities are there- 
fore not easily extrapolated to higher frequencies by adopting 
a spectral index. 

S uch an absolutely calibrated survey was proposed by 
iRei ch & Wielebinski (2000) and carried out at the Dominion 
Radio Astrophysical Observatory (DRAO) in Penticton, 
Canad a, during two observing runs (Ph.D. thesis: IWollebenI 
I2OO5I) . In this paper we give the polarization data, absolutely 
calibrated relative to the Dwingeloo scale and attached to the 
Effelsberg main-beam temperature scale. This northern sky 
survey will ultimately be combined with a southern sky po- 
larization survey a t the same frequency now completed by 
iTestorietalJlEoa. 



2. Receiving system 

The new polarization survey was carried out using the 25.6 m 
telescope at DRAO. The accuracy of the antenna surface makes 
observations at frequencies from about 400 MHz to about 
8.4 GHz possible. At 1.4 GHz the aperture efliciency is 55% 
(10 Jy/K). In the configuration used for this survey the tele- 
scope is equipped with an uncooled L-band receiver operating 
in the frequency range from 1.3 GHz to 1.7 GHz. The tele- 
scope has an equatorial mount with the receiver placed at its 
prim ary focus on three support-struts. The pointing accuracy is 
< 1' llKneell997l) . Various telescop e parameters quoted in th is 
paper were carefully determined bv lHiggs & Tapping! ( I2OOOI) . 
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Figure 1. Block diagram of the 1.4 GHz receiver and contin- 
uum backend of the DRAO 25.6-m telescope used for this sur- 
vey. The polarimeter consists of two analog multipliers as indi- 
cated by cross symbols. 



2.1. Front-end 

The L-band receiver in its original design consists of a 
corrugated-flange feed (scalar feed), which is a scaled copy 

' A fully sampled total intensity survey of the north- 
ern sky at 1.42 GHz with angular resolution of 35' and 
about 15 mK rms-noise. Access to the data is possible via: 
^http://www.mpifr-bonn.mpg.de/survey.html 



of an Effelsberg design JWohlleben etalJll972h followed by 
a piece of circular waveguide and a circular-to-square transi- 
tion. A dual polarization coupler (polarizer) splits the initial 
input signals into its linear (X, Y) polarization components: 
Ey: - aie"^' and = a2e'^'^''^^\ with real amplitudes oi and 02, 
the phase difference 6 between the two components, and ait the 
product of circular frequency and time. 
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Cross-correlation of linear polarization components does 
not yield full measures for the polarization state of a linearly 
polarized incoming wave. Therefore, a quadrature hybrid was 
incorporated to transform linear into circular (R, L) polariza- 
tion components, prior to correlation. The hybrid was placed 
between the feed and the low-noise amplifiers (LNA). This en- 
sures that gain and phase fluctuations of the amplifier chain 
affect only circular hands of polarization, which allows more 
accurate measurement and calibration of Stokes U and Q. The 
outputs of such a quadrature hybrid are phase-shifted sums of 
the two inputs: 

El = 0.5 (flie'("'-f) + a2e"-'^'^^^) , 
with the correlation of both yielding: 

RL = Re(ErE*) = i 0102 - -f sinlip 
— > Stokes U 

LR = Re (Er,-!iEf) = 5 (flj - aj) = f cos 2ip 
— » Stokes Q, 

in case of initially linearly polarized waves, for which the rel- 
ative phase differences is 5 = 0° and the two components are 
a I - aocoscp and 02 = aosintp with the polarization angle 
If and the amplitude aq of the initial wave. Here, RL is the real 
part of the product of the two hands of polarization with the star 
indicating complex conjugate. LR is obtained by phase shifting 
R by 90°, which is done within the polarimeter. For an equato- 
rially mounted telescope, like the DRAO 26-m, the parallactic 
angle of the feed does not change unless the feed is rotated 
with respect to the telescope. For this survey the angle of the 
feed relative to the telescope was fixed at an arbitrary angle, 
and the cross-correlation products required rotation by a con- 
stant angle to be transformed into U and Q in the equatorial 
co-ordinate system. 

2.2. Polarimeter 

The IF-polarimeter is an analog two-channel multiplier It was 
brought from the Max-Planck-Institut fiir Radioastronomie in 
Bonn and is of the type used on the Effelsberg 100-m tele- 
scope for narrow band polarimetry at 1.4 GHz and 2.7 GHz 
jReich et aklliosi) . It operates at an IF of 150 MHz and al- 
lows a maximum bandwidth of 80 MHz, of which 12 MHz are 
used here. The polarimeter provides four output channels: two 
total power channels (RR, LL), and two cross-products (RL, 
LR). To handle quadratic terms and DC offsets caused by the 
analog multipliers and affecting the two cross-products, an in- 
ternal phase shifter is switched with a period of 4 s between 
zero and 180°. Phase shifting one of the two hands of polariza- 
tion (either E,- or Ei) inverts the sign of the correlation products 
but not of the quadratic error term. Hence, correct products are 
given by the difference of the output levels at the two phases: 
RL = RL180 - RLo and LR = LRigo - LRo. 



2.3. System set-up 

Adjustment of the relative gain and phase of the two linear 
hands of polarization ahead of the quadrature hybrid is nec- 
essary. This can be seen by introducing the complex gains 
g i - G^e'"' and gy = Gye'^y with the power gains G.i and G,., 
and the phase mismatch Asxy - Bx- Sy, so that the hybrid out- 
puts read: 

E,- = 5 (fli Gx e'("'+^') + fl2 Gy e'('^'+fv-^/2)^ 
E, = i (fl 1 Gx e'('^'+^' + a2 Gy e'('-'+^v)) 

The effect of mismatches becomes visible in the cross- 
products: 

RL oc a^GxGy sin 2(f cos Ae^v 

LR oc al [gI - G] + [gI + G^) cos 2(p) . ^'^^ 

Hence, phase and gain mismatches of the linear polarization 
components ahead of the hybrid result in non-circular response. 
A similar calculation shows that such mismatches in the circu- 
lar hands of polarization, behind the hybrid, do not affect the 
sensitivity and are corrected by calibration. 

The complete receiving system used for the survey obser- 
vations is shown in Fig.^ The polarizer is followed by direc- 
tional couplers through which a calibration signal is injected 
into the X- and Y-lines. These two lines are cross-coupled by 
the hybrid. Circulators prevent backscattering and reduce in- 
strumental polarization. The R and L hands of polarization are 
amplified by two LNAs^, one for each hand. The signals are 
subsequently filtered, attenuated, and amplified. The IF is cen- 
tred at 150 MHz. The band-limited IF signals are fed into the 
polarimeter backend, which performs square-law detection and 
phase-shifted multiplication of the individual R and L signals. 
The four polarimeter outputs are read by the data acquisition 
PC holding an interface card. 

The calibration signal is generated by a noise source^. The 
period and duration of "cal"-signal injection is different for the 
first and second observing run of this survey. For data obtained 
during the first run the cal was injected every 24 s for a duration 
of 400 ms. For the second run the intensity was reduced and the 
period and duration were changed to 4 s, and 4 s, respectively. 
Also some cables were replaced, which changed the relative 
phase between the two polarization components of the cal and 
thus its polarization angle. 

The system temperature (Tsys) of the receiver without sky 
emission was measured from ground radiation profiles and 
amounts to 125 + 10 K in total power at the zenith. The fol- 
lowing main contributors to T^y^ are: 1) thermal noise of the 
LNAs contributes ~ 35 K; 2) the hybrid and additional coax 
cables have a loss of about 0.45 dB and raise Tsys by ~ 30 K; 3) 
the circulators have an insertion loss of 0.25 dB and thus add 
~ 17 K; 4) the directional couplers have an insertion loss in 
their main line of 0.2 dB which adds 15 K; and 5) additional 
15 K allow for ground radiation, atmosphere, and other noise 
sources. 

- Berkshire L- 1 .4-45HR with quoted gains of 57.0 dB, and 5 1 .6 dB, 
respectively 
3 Noise/Com Inc Model NC3 10 IE 
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Table 1. Receiver and antenna specifications. 



Telescope coordinates 


-119° 37.2' 




+49° 19.2' 


Antenna diameter 


25.6 m 


HPBW (effective) 


36' 


Aperture efficiency 


55% 


Bandwidtii (November 2002) 


12 MHz (10 MHz) 


Pointing accuracy 


< 1' 


Intermediate frequency 


150 MHz 


System temperature 


125 + lOK 


Gain variations in 30 days 


< 4% 


Piiase tracking across band 


~ 5° 



2.4. Software 

Software for data acquisition, reduction, and calibration was 
written in glish, a programming lan guage that is part of 
the AJ PS++ environment (e.g.: Wini3 l2002l iMcMulhn et alJ 
12004 1). We utilized two PCs for the backend control and data 
acquisition, and an additional IBM 520 computer for antenna 
control and scheduling of observations. The data acquisition 
PC was equipped with an interface card^ carrying four coun- 
ters that read out signal levels at the polarimeter output lines. 
The integration time for a single integration was 40 ms. On 
the data storage PC several data streams carrying packed inte- 
grations, information about the current telescope status, actual 
observations, telescope position, local sidereal time (LST), and 
wind speed are merged on arrival and stored onto disk. 

3. Observing strategy 

Observations were made in two sessions, from November 2002 
through May 2003 (first run) and from June 2004 through 
March 2005 (second run). The central frequency was set to 
1.41 GHz, avoiding the Galactic Hi emission. Observations 
covered the declination range +90° to -29°. The location of 
drift scans is shown in Fig. |2] To minimize the influence of 
varying ground radiation, all observations were made as drift 
scans with the telescope stationary on the meridian. Drift scans 
were scheduled as subscans starting every 60 min, with 57 min 
length each, to allow position control of the antenna between 
subscans. The observations were carried out with the telescope 
left unattended, whereas elevation scans and drift scans at de- 
clinations > 85° were done manually. Declinations were ob- 
served in random order to reduce systematic effects. Drift scans 
were started (stopped) 1 to 2 hours before sunset (after sunrise). 
Each drift scan was fully sampled according to the Nyquist cri- 
terion (0.5 X HPBW). Over the declination range -29° to 90° 
we achieved 41.7% of full Nyquist sampling in the telescope 
time allocated for this survey. 

3.1. NCP observations 

The circularity of the system response to polarized emission 
was checked using measurements of the Northern Celestial 
Pole (NCP), because the rotation of the polarization angle rel- 

National Instruments: NI 6601 
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Figure 2. Coverage of the sky observable from DRAO with sur- 
vey drift scans (black lines). Total coverage of the region shown 
is 41.7% of fufl Nyquist (0.5 x HPBW) sampling. This figure 
is also available from the web-page (see Sec.|6} and the exact 
declinations of the drift scans can be taken from that. The gaps 
apparent at some declinations and right ascension intervals are 
due to the incomplete coverage and a randomized observing 
strategy. 



ative to the stationary telescope can be measured at that point. 
A calibrated 1.5 day observation of the NCP is displayed in 
Figure|3l Based on this scan the polarized intensity of the NCP 
is 80 mK, which is, within the errors, in agreement with the 
temperature of 60 mK at 1.4 GHz derived by the LDPS. The 
polarization angle is -37° in RA-DEC coordinates. 
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Figures. A 1.5 day observation of the NCP, obtained on May 14th through May 16th, 2003. Shown are measurements made at 
daytime (20 h - 1 1 h LST) and by night. The soHd Hne shows a sine wave fitted to the data, taking a Hnear basehne into account. 
These data were caHbrated in the following way: the electronic gain was corrected using the noise source (Sect. lTm and the 
Mueller matrix (Sect. l433t was applied. 



3.2. Ground radiation effects 

At centimeter wavelengths rough and dry ground behaves like 
a black body of appr oximately 250 K effective temperature 
llAnderson et alJlTQQll) . If the surface is smooth compared to 
the observing wavelength it also reflects emission from the sky. 
Moreover, ground radiation off'sets in the data can be expected 
to be time variable. Seasonal as well as day-night variations 
of the ambient temperature (Tamb) lead to fluctuations in the 
system temperature. In addition, the level of ground-reflected 
signals picked up by the sidelobes depend on LST while drift 
scanning the sky. 

In the following we estimate the amplitude of the summer- 
winter variation of the ground radiation offsets in total inten- 
sity. If the antenna is pointed at an elevation of ~ 20°, half 
of its sidelobes receive radiation from the ground while the 
other half and the main lobe see the atmosphere. With the ra- 
diative efficiency t/r = 0.995, the stray factor /3 = 0.26 (the 
fraction of the antenna solid angle contributed by responses 
outsi de the main beam), and the spillover efficiency rj^po = 
0.95 llHiggs & Tappin J200nl) . the sidelobe contribution to the 
ground offsets in total power is rjjtfiiT^^b + T^itn)/! and the con- 
tribution of spillover noise is ?7r(1 - 77spo)(7^amb + 7^atm)/2. For 
the vertical telescope (~ 90° elevation) these are (0.1 Tamb + 
0.9 ratm)?7Ry6, and 77r(1 - /7spo)ramb, assuming that the back 
lobes contribute 10% to the total signal r eceived through the 
sidelobes. With T^^u = 2 K (Gibbins 1986) and Ta^b = 250 K 
and 240 K during summer and winter, respectively, the seasonal 
variation of the ground radiation in total power is 1 .6 K for the 
horizontal, and 0.7 K for the vertical telescope. 

Ground radiation is received through the side and back 
lobes of the antenna and raises the system temperature, which 



is observed as variations of the baselevel (offsets). Numerical 
calculations show that the far sidelobes are polarized by up to 
50% (INg et alJl2005l) . Therefore, ground radiation may cause 
offsets in the cross-products as well as in total power, with am- 
plitudes of several hundreds of mK. On the basis of the mea- 
surements made for this survey it is impossible to judge how 
much the ground radiation itself is polarized, or whether the 
relevant sidelobes are polarized. 

4. Data processing and calibration 

In the following we summarize the data processing steps ap- 
plied to the raw data. A flow chart of the reduction chain is 
displayed in Fig.|3 Here, we neglect the total power data (RR, 
LL), which are only used for the correction of main beam in- 
strumental polarization during data reduction and are not in- 
tended to be published. 

4.1. General overall Strategy 

For the absolute calibration of the Stokes U and Q data the fol- 
lowing calibration scheme was applied. The LDPS at 1.4 GHz 
provided 1666 absolutely calibrated pointings observed north 
of 0° declination. We used a digitized version of these data, 
kindly provided by T. Spoelstra, as reference values for the de- 
termination of the system parameters of our instrument and the 
determination of zero levels in U and Q. While surveying the 
sky, 946 of these pointings (referred to as congruent pointings if 
their centers match within a radius of 15') were observed in the 
course of the drift scans and used for the caUbration. In this way 
we efi'ectively corrected our observations north of O°declination 
for ground radiation. For the extrapolation of the zero levels be- 
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Figure 4. The data reduction chain as explained in the text. 



low 0° declination, profiles of ground radiation were measured. 
In a final calibration step, the main-beam brightness tempera- 
ture scale in U and Q of the DRAO survey was refined by com- 
parison with data taken from the EMLS (see Section !?. 3. 7> . 

Serious time-variable errors were found to arise from fluc- 
tuations of Tsys- These could be corrected by comparing drift 
scans with neighboring scans and thus separating instrumental 
variations from variations of the sky brightness temperature in 
U and Q. Fluctuations of the electronic gain, most likely caused 
by the LNAs, were corrected using the cal. Solar interference 
and ionospheric Faraday rotation could largely be avoided by 
not observing during day-time. Other time-variable errors such 
as radio-frequency interference (RFI) or weather changes were 
recognized and flagged during data reduction. 

4.2. RFI flagging 

The remote DRAO site is surrounded by mountains which 
block RFI to some extent. In addition, most of the institute 
buildings are screened by a Faraday cage, which reduces most 
of the RFI caused by local facilities. Therefore, DRAO pro- 
vides ideal conditions for radio observations in L-band. In fact, 
only 0.1% of the total observing time was lost due to inter- 
ference. Most of the RFI observed was of short duration and 
distorted only a few integrations. 

Interference was detected by an algorithm that searches for 
short-term peaks separately in all four channels. This search 
was performed within a 6 minute time window, which was 
moved along the time-axis of the raw data in steps of 90 s. 
Integrations with signal levels higher than 8 o-i-ms (total inten- 
sity) or 5 (Tj^s (polarization) above the average signal level in 
one or more channels within the window were flagged. This 
RFI flagging was done a second time on averaged data with 
a threshold of 3 cTi-ms- A third, visual flagging was done after 
gridding of the data. All four channels of flagged integrations 
were removed from further data processing. 



4.3. Calibration and gridding 

The following steps perform a calibration of both the relative 
intensity scale and the absolute offsets in Stokes U and Q. 



4.3.1 . Correction for electronic gain 

The calibration noise source provides a constant polarized sig- 
nal whose actual intensity and angle are extracted from the 
raw data by separating the amplitudes in the auto- and cross- 
products from sky emission. Any change of Plcai or PA^^i is 
interpreted as due to gain or phase fluctuations within the re- 
ceiver The strength of the cal is estimated to be (Plcai) = 36 K 
and (PAcai) = 71° during the first observing run, and (Plcai) - 
1 K and (PAcai) - -50° during the second run, where the angle 
is given with respect to the NCP, measured from north through 
east. 

The raw data are corrected for electronic gain by referring 
the recorded intensities of the four channels to the cal values: 



(5) 



RR' = RR/<RRcai) 

LL' = LL/(LLeal) 

RL' = (RL cos 2 (PAcai) -LR sin 2 (PAcai)) /(Plcai) 
LR' = (LRcos2(PAcai) + RLsin2(PAcai))/(PIcai), 

with (Pleal) = ((RLeal)' + (LRcal)')'^' and (PAeal) = 0.5 X 

arctan((RLcai) / (LRcai)). For data obtained during the first ob- 
serving run, integrations that include the cal were flagged after 
gain correction because of the intensity of the cal, which re- 
duced the effective observing time by 1.7%. 

4.3.2. Determination of tine IVIueller matrix 

To correct (time-invariable) systematic errors, the raw data 
were calibrate d via the inverse of a system Mueller matrix 
(Mu elleJl948l) . which was derived in two iterative steps using 
the LDPS as reference. The 4x4 Mueller matrix of a receiving 
system describes most generally the transformation of the in- 
coming signal power by the receiver components: Sm,, = Al-5in 
where S is any of the Stokes vectors (see e.g.: iHeiles et alJ 
Heiles 2Q Q2i|). Here, the inverse Ai' was derived by fit- 
ting its entries to the equation formed of the 4-vectors of the 
(electronic-gain-corrected) DRAO raw data ("raw") and the 
reference values ("ref"): 



(6) 



Differing from the conventional definition, the matrix trans- 
forms correlation products instead of Stokes parameters, which 
are converted according to: RR = I + V, LL = / - V, RL - U, 
and LR - Q. The conversions RL - U and LR - Q are possi- 
ble because the parallactic angle of the feed is constant. 

The intention of this first Mueller-matrix correction is to 
convert the arbitrary units of the polarimeter outputs (counts) 
into a rough temperature scale. The total power channels (RR, 
LL) are calibrated to approximate brightness temperature us- 
ing the Stockert 1 .4 GHz total intensity (Stokes /stockert) survey. 
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The cross-products (RL, LR) are scaled and rotated to roughly 
express Stokes U and Q measured in Kelvin, using the LDPS. 
This represents the first iteration step of a two-step iterative de- 
termination of the inverted instrumental Mueller matrix. 

The reference values for the computation of the pre- 
calibration matrix are given by: 



RRrf = (0.5 X 1.55 (/stockert - 2.8 K)) + 2.8 K 

LL,ef = (0.5 X 1.55 (/stockert - 2.8 K)) + 2.8 K 

RLref = t^LDPS 

LRref - 2lDPS> 



(7) 



which includes the transformation of full-beam into main-beam 
brightness temperature of the Stockert survey with a scaling 
factor of 1 .55 and t he correction for extr agalactic background 
emission of 2.8 K jReich & Reichlll988h . which is subtracted 
and subsequently added. Stokes V is assumed zero. 

The fitting of the pre-calibration matrix proceeds using 
Eq.|6l For each point, Stockert data supply RRief - LLref. The 
Dwingeloo telescope provides RL,ef and LRief . The DRAO 26- 
m telescope provides all matrix elements on the right hand side 
of Eq. |6l allowing M\ to be determined. Hence, if the total 
number of pointings is A^, the resulting set of equations with 
n = l,2,...,AforRRandLLreads: 



RR,ef,„ = miifRR;ky.„+RR;ff,„) 

LLref,« = »I22 (LL^j^^ „ + LL^jf . 

For RL and LR it is: 

RUef,„ = m33(RL;k,,„+RL;J + m34(LR:,^,„ + LR;^„) 

LRref,„ = m43(RL^ky.„+RUff,„)+'«44(LR:k,,„+LR;^„). 



(8) 



(9) 



At this point the raw data still contain time-variable offsets 
("off") in addition to the sky values ("sky"). Therefore, the ma- 
trix entries oth, mii, m33, 11144, '«34 and m43 of the matrix M[ 
can only be fitted iteratively by approximation of these offsets 
as done in the following algorithm: 

1 . The DRAO raw data for congruent pointings are corrected 
for electronic gain. 

2. The initial matrix Ai[ is set to be the unit matrix. 

3. The DRAO raw data are corrected by multiplication with 
M[ according to Eq.|6] 

4. Ground radiation and system temperature drifts are ap- 
proximately removed by subtracting a linear baseline fitted 
through all congruent pointings at equal declination from 
0*" to 24*^ right ascension. 

5. The fit-error is given as the sum of the squared differences 
between corrected DRAO data and the reference values. 

6. Matrix entries of M[ are either randomly altered and the 
iteration continues at point three or the iteration is stopped 
if the minimization of the fit-error is achieved. 

This way, pre-calibration matrices have been determined for 
each observing month, using about 100 congruent pointings 
per month. These monthly matrices were checked for time- 
variability and then averaged. No systematic variation of the 
matrix entries with time were found bigger than the fitting er- 
rors. 



4.3.3. Application of tine Mueller matrix 

We used the following matrix for the pre-calibration of raw data 
obtained in November 2002: 



M' 



K/ counts 



17 o^ 

96 
-18.0 43.7 
56.0 17.2 



(10) 



with errors of roughly 20% in the matrix entries. Two other ma- 
trices were used to describe the slightly modified receiver dur- 
ing periods from December 2002 through May 2003 (change 
of IF filters in November 2002) and June 2004 through March 
2005 (replacement of hybrid and adjustment of cal signal in 
May 2004). These matrices reveal changes of the system gain 
and phase, caused by the modifications mentioned above. Raw 
data were thus pre-calibrated according to Eq.|6lusing the ma- 
trix corresponding to the time of observation. 

4.3.4. Correction for ground radiation 

Profiles of ground radiation were measured by sweeping 
the telescope between -30° and -1-90° declination along the 
Meridian. A series of elevation scans were made at different 
LST to average out emission from the sky. We used four such 
scans, made during the night of 2002 December, 12th/13th. A 
complete scan took 1.5 h with a scanning speed of 1.3°/min. 
Elevation scans were corrected for electronic gain fluctuations 
and converted into brightness temperature by applying the pre- 
calibration matrix. Ground radiation profiles for each channel 
are determined on the basis of these elevation scans. 

Figure Is] shows the derived ground radiation profiles. The 
mean offsets of about 3.5 K in the Stokes U profile and -3.7 K 
in the Q profile may indicate: 

1 . remaining phase and gain mismatches in the receiver pre- 
ceding the hybrid, or 

2. correlated receiver noise caused by cross-talk, or 

3. polarized backlobes of the telescope, picking up ground ra- 
diation. 

The first can be excluded from system checks; more likely are 
the second and third. Also interesting is the north-south asym- 
metry about 90° elevation of the U and Q profiles, which is 
probably due to the multipole structure of the response pattern 
of the telescope to unpolarized radiation. The response pattern 
of the DRAO 25.6 -m tele scope in total power was measured by 
Higgs&Tanninll ll200nl) and shows sidelobe structure to radii 
as far as 45° off the main beam at a level above -49 dB. These 
are mostly caused by aperture blockage due to the three support 
struts and the receiver box. Irregularities of the antenna surface 
complicate the sidelobe structure. 

Subtraction of ground radiation from the raw data is 
straightforward. The ground offsets for each declination (be- 
tween -29° and -1-90° declination) and channel are taken from 
the profiles. These offsets are then subtracted from the drift 
scans. This way, the declination-dependent component is re- 
moved and only the time-variable components of ground radi- 
ation and system temperature are left; these are corrected later 
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Figures. Ground radiation profiles for Stokes U and Q determined on the basis of elevation scans. A possible explanation for 
the north-south asymmetry is given in the text. 



(see Section l4.3.6> . The remaining uncertainties of this proce- 
dure are around 50 mK or less, according to the variations be- 
tween individual elevation scans at positions where polarized 
emission from the sky is low. 

4.3.5. Solar interference and ionospheric Faraday 
rotation 

The data were visually inspected by comparing drift scans 
with neighboring scans to identify high levels of ionospheric 
Faraday rotation or solar interference. Only some of the day- 
time data, observed shortly before sunset and after sunrise, 
show apparent solar interference and needed to be flagged. No 
indication was found for ionospheric Faraday rotation at night. 

4.3.6. System temperature fluctuations 

Data reduction was made more difficult by fluctuations of sys- 
tem temperature, Tsys. Fluctuations of the order of several parts 
in 1000 (a few hundred mK in ~ 100 K) were subsequently 
found to have occurred throughout both observing periods. 
Fluctuations in the RHCP and LHCP channels were not corre- 
lated, and there was no apparent correlation of the fluctuations 
of Jsys with time of year, temperature, rain, or other variables. 
Over periods of months we found the amplitude of these vari- 
ations to be several Kelvin. In subsequent observing nights the 
baselevels of drift scans dififered by several hundreds of mK 
(after correction for ground radiation). Such fluctuations oc- 
curred on time scales from hours to months. The origin of these 
fluctuations is unknown. 

The Tsys-fluctuations must be corrected while large-scale 
sky emission in the data should be preserved to obtain abso- 
lutely calibrated maps. However, variations of the system tem- 
perature cannot easily be distinguished from real sky emission. 
An algorithm was applied that effectively separates random 
fluctuations from systematic sky emission (see App. |AJ. The 



remaining errors are of the order of 50 mK at maximum as in- 
dicated by simulations^. 

4.3.7. Refinement of the Mueller matrix 

A more accurate scaling and rotation of the RL and LR scales 
and a correction for main beam instrumental polarization were 
achieved by the second calibration. No further corrections of 
the total power data were made as the intention was the cali- 
bration of polarization. Instead of iterative fitting, here least- 
square fitting can be applied because the raw data are now 
corrected for ground radiation offsets. The reference points are 
again taken from the LDPS. 

A single calibration matrix is determined for the whole data 
set. The following correction matrix was derived: 

(I 0^ 

M' = 10 

2nd 0.001 0.0001 1.125 -0.013 ' ^ ' 

.-0.009 -0.011 0.051 1.240. 

with errors of around +0.01 in the matrix entries. Raw data 
were corrected according to Eq. |6] The second calibration 
matrix mainly corrected scaling errors that were introduced 
by the separation of Tsys-fluctuations. Instrumental polariza- 
tion within the main-beam (cross-talk) is given by RLinst - 

mini.'i 1 • RR + OT2nd,32 ' LL and LRi„st = 'W2nd,41 • RR + 'W2nd,42 " LL 

and hence amounts to 0.1% in Stokes U and -2.0% in Q. 

4.3.8. Gridding and interpolation 

Integrations were averaged and binned into a grid of equatorial 
coordinates with 15' cell size, which is less then the telescope's 
HPBW/2 required for full Nyquist sampling. No weighting 
scheme was used; measured values were uniformly weighted 
and averaged across the cell. 

^ For this, the algorithm was applied to artificial Stokes U and Q 
maps with simulated noise and Tgyj -fluctuations added to them. 
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An interpolation between scans needs to be applied. In con- 
trast to "uniform" undersampling, the DRAO survey is fully 
sampled along right ascension and undersampled in declina- 
tion. Considering this specific sampling problem an interpo- 
lation routine working in the Fourier domain was applied. 
The drifts scans were Fourier transformed, using the DFT^- 
algorithm of AIPS++. Unobserved declinations were then filled 
by linear interpolation in the Fourier space. The interpolated 
scans were then reverse transformed into the image plane. 



Table 2. Survey specifications. 



Integration time per 15' in R.A. 
Observing period 

Final rms-noise 
Systematic errors 
Declination range 
Fully sampled area 



60s 

November 2002 - May 2003 

June 2004 - March 2005 

12 mK 

< 50 mK 

-29° to +90° 

41.7% 



4.3.9. Temperature scale refinement with tine ElVILS 

Brightness temperatures from the DRAO survey were com- 
pared with the EMLS over an area of 1800 square degrees 
of sky. Two preliminary maps from the EMLS were used: a 
35° X 35° map towards Cassiopeia-A and a 40° x 16° map cov- 
ering the Taurus-Auriga-Perseus molecular clouds. 

As the EMLS lacks large-scale emission a direct compari- 
son with the DRAO survey is not meaningful. Therefore, miss- 
ing large-scale structures were approximated by polynomials. 
The coefficients of these polynomials and the scaling factor of 
the DRAO temperature scale were fitted to the Effelsberg data 
over the entire region, with the difference between EMLS and 
DRAO survey to be minimized. Figure |6l shows four constant 
declination scans, two through each comparison region. 

The calibrati on of the EMLS is based on standard cali- 
bration sources jUvaniker et akl 'lQQS). particularly the well- 
established value for the polarized flux density of 3C 286, and 
can thus be assumed to be more accurate than the LDPS cal- 
ibration. The LDPS is not sufficiently sensitive to allow cali- 
bration by compact sources, but relies on a set of calibration 
points (a = 57°0, 6 = 64?0 and a = 240?0, 6 = 23°0) on 
the sky where polarized intensit ies were d etermined earlier at 
Cambridge with a larger beam ( Binghaml[l966 ). Frequent ob- 
servations of these calibration points gave a high internal con- 
sistency to the LDPS, but the accuracy of the scale is not com- 
parable with that achievable at the Effelsberg telescope today. 

5. Accuracy and errors 

The theoretical rms noise for the correlation receiver used here, 
with a bandwidth of 12 MHz, and an integration time of 60 s, 
calculates to 3 mK (valid for Stokes U and Q). In the follow- 
ing we summarize the measured rms noise and the systematic 
errors of this survey. 

5.1. RMS estimation 

Taking the NCP measurements (Sect. 13. 1> . the deviation from 
a circle in the U-Q plane gives the resulting rms noise of the 
DRAO survey. This is, for an integration time of 60 s, 12 mK 
rms in Stokes U and Q taking only night-time data into account, 
and ~ 35 mK including day-time data. 

Figure shows the correlation of the final data with the 
reference values, prior to adjusting the DRAO intensities to 
the EMLS. The correlation coefficients are ru - 0.89 and 

^ Discrete Fourier Transformation 



rg = 0.86. Numerical calculations^ show that, to reproduce 
these correlation coefficients, the rms noise in the DRAO data 
must be 12 mK in Stokes U, and 33 mK in Q. The differ- 
ence may reflect an error in the Dwingeloo data higher than 
the quoted error of 60 mK. We note that, based on these cal- 
culations and the LDPS data, the best achievable correlation 
between the two data sets is ru - rq - 0.90. 

5.2. Resulting error 

Based on the NCP measurement and the correlation coefficients 
the final rms noise in Stokes U and 2 is 12 mK, which is 
9 mK higher than the theoretical rms noise. Systematic errors 
are introduced during the separation of sky emission and T^ys- 
fluctuations. In case of intense polarized structures, which are 
elongated along right ascension, this error may be < 50 mK, de- 
pending on the shape, intensity, and coverage of the structure. 
At decUnations < 0° a systematic baselevel error of < 50 mK 
may be caused by uncertainties in the ground profiles, which 
should be highest at the lowest declination. 

5.3. Instrumental artifacts 

A correction for instrumental polarization caused by sidelobes 
is not intended for this survey. This requires fully sampled total 
intensity maps and precise measurements of the response pat- 
tern of the DRAO telescope in polarization. Some instrumental 
polarization effects can be seen around intense compact sources 
in Figures El and |9l (e.g. Taurus-A at RA = 5*' 35™, DEC = 
-5° 30' in Figure|5). Instrumental effects are also seen around 
the intense emission from the Galactic plane in directions to- 
wards the inner Galaxy (from RA = 11^ 38™, DEC = -29° to 
RA = 19*' 40™, DEC = 20°) seen especially in the U image 
of FigurelH These instrumental effects do not extend beyond a 
few degrees. Some instrumental effects have the appearance of 
scanning artifacts, clearly apparent in the image in Galactic co- 
ordinates around the North Celestial Pole at / - 123°, b - 27°. 
These are due to undersampling in this region. 

Based on a number of point sources we estimate the instru- 
mental polarization of the first sidelobe to be < 6% at maxi- 
mum. For diffuse, extended emission the effect of sidelobes in 
most likely negligible because their contributions over a large 
area around the main beam average close to zero. 



' These assume the quoted error of 60 mK of the LDPS, and a nor- 
mal distribution of U and Q around zero with a FWUM of 120 mK. 
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Figure 6. Four scans through the two test regions used to establish the calibration of the survey relative to EMLS observations 
(thin solid lines), two from the Taurus complex (left panels) and two from the Cassiopeia-A region (right panels). DRAO survey 
data are drawn as thick solid lines. The large-scale emission missing from the EMLS data are shown as polynomials (whose 
characteristics were determined by an iterative technique and added to the EMLS data shown in this figure - see text) and are 
drawn as dashed lines. The declinations of the scans are indicated. 



6. Final maps and data availability 

The final maps in Stokes U and Q are displayed in Figures |8l 
and|9lin equatorial co-ordinates in the rectangular projection in 
which the data processing and gridding were done. Figure [Tol 
shows polarized intensity in Galactic co-ordinates in Aitoff pro- 
jection. In this section we will comment on the most prominent 
astronomical results. Detailed interpretation of the newly dis- 
covered features must await a forthcoming paper. 

The most intensely polarized features of the Northern sky 
are the North Polar Spur, at Galactic latitudes higher than 25°, 
extending from longitude 300° to 40°, and the Fan Region, of 
extent 60° by 30°, centred on I = 140°, b ^ 5°. These features 
were readily apparent in older data (e.g. iBrouw & Spoelstr^ 
Some large nearby H ii regions can be seen as depolar- 
ization zones. Particularly obvious are Sharpless 27 (/ = 6°, 
b = 23°), distance 170 pc, and Sharpless 220 (/ = 163°, 
b = -15°), distance 400 pc. Such depolarization effects can be 
used to establish a minimum distance to the polarized emission 
regions. 



The North Polar Spur coincides with a bright emission fea- 
ture seen in total intensity in many Galactic surveys. It is gen- 
erally considered to be a supernova remnant (SNR), generated 
by the explosion of a mem ber of the Sco-Cen association at 
a distance of about 150 pc / Berkhuii serJI 1 97 ll) . In Fig.llllwe 
show the southern part of this spur. 

The Fan Region, in contrast, has no obvious counterpart in 
total intensity. A polarization map is shown in Fig. ^2 Earlier 
investigators considered this polarized emission to originat e at 
distances under 500 pc ( Wilkinson & Smith 1974; Spoelstral 
Il984h . Hn regions seen in the same direction did not appear 
to depolarize the emission from the Fan Region, and the po- 
larized emission was therefore considered to be closer than the 
H n regions, whose distances were known. With our new data, 
which has greatly improved sampling and sensitivity, we can 
re-examine this conclusion. We have detected definite depolar- 
ization by a nuni ber of H ii regions seen against the Fan Region 
lWolleben"2005"). On the basis of this new evidence we con- 
clude that the polarized emission from the Fan Region orig- 
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Figure 7. Correlation of Stokes U (left panel) and Q (right panel) values from the DRAO survey with the Dwingeloo survey 
before adjusting the temperatures to the EflFelsberg scale. 



inates over a range of distances, extending from 500 pc to a 
few kpc, the latter distance corresponding to the distance of the 
Perseus arm. Over much of the Fan Region the electric vectors 
are perpendicular to the Galactic plane, indicating that the mag- 
netic field is aligned with the plane. The fractional polarization 
in this region is up to 50%, not much lower than the theoretical 
maximum of 70%, indicating that there must be a very regular 
magnetic field structure and very little depolarization along a 
long line of sight, extending right through the Perseus arm. 

A new result that is apparent in our data is a depolariza- 
tion region that extends between latitudes +30° and -30°, from 
longitude about 65° towards longitude zero. Our superior cov- 
erage as well as sensitivity and the representation of the data 
in grayscale, as in Figure^] reveals the small-scale structure 
in this region. The typical feature size in this region is no more 
than a few degrees. Earlier surveys with their sampling of order 
two degrees could not show this region adequately. The depo- 
larization must arise in a very nearby magneto-ionic region: the 
emission from the North Polar Spur is depolarized, and the dis- 
tance to the Spur is no more than 150 pc. Particularly striking 
are the sharp upper and lower boundaries of this depolarization 
region. 

Examination of polarization angles at Galactic latitudes 
above 70°, both north and south, shows systematic changes 
with position and a remarkable symmetry of polarized inten- 
sity. Fractional polarization reaches values of 50% (northern 
pole) and 30% (southern pole). These results suggest that the 
Sun is located inside a synchrotron emitting shell with the mag- 
netic field parallel to the Galactic plane. 

The 1.41 GHz polarization data presented here are avail- 
able in Galactic or equatorial coordinates via the following 
web-pages; http://www.mpifr-bonn.mpg.de/div/konti/26msurvey 
or http://www.drao.nrc.ca/26msurvey. Alternatively, the 



reader may contact MW at maik.wolleben@nrc-cnrc.gc.ca. 
Ungridded data may be made available on request in form of 
tables. 



7. Summary 

We have presented a new survey of linear polarization of the 
northern sky, obtained using the DRAO 25.6-m telescope at 
1 .4 GHz. Observations were made by drift scanning the sky to 
keep ground radiation as constant as possible, which allows an 
absolute calibration of the Stokes U and Q maps. A fully sam- 
pled area of 41.7% has been observed. Although a complete 
coverage of the northern sky could not be achieved, this data 
base provides about 200 times more data points than data so 
far available and sensitivity 5 times better than previous data. 
The new survey reveals previously unknown structures and ob- 
jects in polarization and provides valuable information for the 
derivation of foreground templates of polarized emission re- 
quired for CMB experiments. 
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Figure 9. Stokes U (top) and Q (bottom) in the right ascension interval from O'' to 12*' in equatorial coordinates. 
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Figure 10. Polarized intensity in Galactic coordinates. The projection center is at / = = 0°. The co-ordinate grid is in steps of 
60° in Galactic longitude and 30° in Galactic latitude. 
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Figure 11. The two maps of polarized intensity, with E-vectors of the polarized emission overlayed, show in greater detail part of 
the Fan-region (left panel) and the southern extension of the North-Polar Spur (right panel). In both maps depolarization caused 
by local Hn-regions is apparent (NGC 1499: from / = 145° to 170°, h < -10°; and Sharpless 27: roughly 10° in diameter, 
centered at / = 6° and b = 23°). The synchrotron emission of the North-Polar spur seems to be depolarized atb< 25° (from 
/ = 30° to 40°). In both maps contours of polarized intensity run from mK to 800 mK in steps of 100 mK. Vectors are shown 
on a grid of 1.25°. 



Appendix A: Algorithm for correction of 
Tsys-fluctuations 

This algorithm compares each drift scan pairwise with its 
neighboring scans. The basic assumption is made that T^ys- 
fluctuations are random on time-scales of weeks, which means 
that fluctuations observed at the same declination but weeks 
apart are assumed to be uncorrected. Therefore neighboring 



drift scans with shghtly different dechnations were observed 
weeks apart so that system temperature fluctuations might dif- 
fer, but not the large-scale structure contained in these scans. 
The declination range from which neighboring scans are se- 
lected is adjustable in the program as described in the follow- 
ing. 

The following performs an iterative separation of random 
and systematic structures in the drift scans. Let m be an index 
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numbering all observed drift scans. The following loop is ap- 
plied until m has reached the total number of scans: 

1. To the drift scan with index m all neighboring scans n 
within +c5° in declination are assorted. Values for 6 between 
3 and 40 were found to be useful here. These numbers re- 
flect the sampling and should decrease with more complete 
coverage or larger number of drift scans, respectively. 

2. Each pair of drift scans (m,n) is convolved with a Gaussian. 
The width of the Gaussian is a times the separation in dec- 
Unation Adm,n = \dm - d„\ of the two drift scans. This re- 
moves spatial structures smaUer than cr x Adm^. Values for 
cr between 4 and 100 are useful here. 

3. The weight Wm_„ - \S - Adm,n\ is introduced. 

4. The difference of each pair of drift scans AT„^„ is cal- 
culated. The sum of these difl'erences multiplied with the 
weight Wm,n is then subtracted from the ;7i-th drift scan: 
T'm = Tm- Y^n ^m.n X Ar„,„ with T denoting brightness 
temperature. 

5. Index m is increased by one and the loop restarted until m 
reaches the number of scans. 

This algorithm was applied several times with different param- 
eters until a satisfactory separation of sky emission and sys- 
tem temperature was achieved. The following parameters were 
used: loops 1 to 4: 5 = 40 , cr = 100; loop 5: (5 = 10 , cr = 15; 
loops 6 - 9: ^ = 6 , cr = 8; loops 10 - 1 1: 5 = 3 , cr = 4. 
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